Casimir force at both non-zero temperature and finite conductivity 
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We find the joint effect of non-zero temperature and finite conductivity onto the Casimir force 
between real metals. Configurations of two parallel plates and a sphere (lens) above a plate are 
considered. Perturbation theory in two parameters (the relative temperature and the relative pen- 
etration depth of zero point oscillations into the metal) is developed. Perturbative results are 
compared with computations. Recent evidence concerning possible existence of large temperature 
corrections at small separations between the real metals is not supported. 
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Currently the Casimir effect is attracting considerable 
interest. A large amount of information is available nowa- 
days of both theoretical and experimental nature. Theo- 
retical progress was made in elaborating different approx- 
imate methods [1-3] and in the problem of a dielectric 
sphere where, for instance, the structure of the ultravi- 
olet divergencies had been clarified Q. In || the addi- 
tive method was successfully applied to a dilute dielectric 
ball, and in Q a progress was made in obtaining analyti- 
cal results. The Casimir force was demonstrated between 
metallic surfaces of a spherical lens (sphere) above a disk 
using a torsion pendulum Q and an atomic force micro- 
scope [p. 

The increased accuracy of Casimir force measurements 
invites a further investigation of different theoretical cor- 
rections. In pi the Casimir force for the configuration 
of a sphere above a plate was computed by taking into 
account surface roughness and finite conductivity correc- 
tions up to fourth order in respective small parameters. 
That result is in excellent agreement with the measured 
Casimir force. Except for contributions of surface rough- 
ness and finite conductivity, corrections due to non-zero 
temperature play a dominant role above some distance 
between the test bodies. The general expression for the 
temperature Casimir force between dielectric plates was 
firstly obtained in Q (see also 12 1). The temperature 
Casimir force between perfectly conducting plates was 
found in [13-15], including the limiting cases of large 
and small plate separations (high and low temperatures). 
These results were modified for the configuration of a 
spherical lens above a disk in 0. The temperature cor- 
rections are found to be insignificant within the separa- 
tions of experiments |||)| (from a «0.1 /im till a=0.9/zm 
or 0.5 fim). As for experiment they constitute up to 
174% of the net force at room temperature T=300K at 
the largest separation a=6 pim [|l6| (in spite of this, exper- 
imental data of Q are not sufficiently accurate to demon- 
strate temperature corrections). In particular, it should 
be emphasized that the joint effect of non-zero temper- 



ature and finite conductivity of the boundary metal was 
not investigated up to the present. 

The computations of the recent paper |l7j have cast 
some doubt on the possibility to use the low temperature 
limit of the temperature corrections for perfect conduc- 
tors [13-15] in order to describe the real metals. For the 
configuration of a sphere above a disk the difference of 
4pN at a ~0.1 /im was found depending on whether one 
uses the expression for the temperature Casimir force 
for a real metal or the zero temperature one. This is 
in contradiction with the generally accepted behavior of 
the temperature correction between perfect conductors 
at small separations which is proportional to (kBTa/Uc) 3 
for a sphere above a disk, ks being the Bolzmann con- 
stant [7,13-15]. The authors of JlTj hypothesized that 
for a real conductor the temperature correction to the 
Casimir force at low temperature can behave as ksTa/hc 
and be important. 

Here we present a perturbative calculation of the joint 
influence of non-zero temperature and finite conductivity 
on the Casimir force. The obtained results arc the gener- 
alization of [13-15] to the case of real metals. They can 
be used for the interpretation of precision experiments 
on Casimir force. No unexpected large temperature con- 
tributions arise at small separations. The erroneous as- 
sumption of p7| is explained below in detail. 

We start with the configuration of two plane parallel 
plates with the dielectric permittivity e separated by an 
empty gap of thickness a. At arbitrary temperature T 
the attractive force per unit area acting between plates 
is given by the Lifshitz formula Jl^ ] 
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The prime on the sum indicates that the term with n = 
is to be taken with the coefficient 1/2. Let us now 
introduce new variables x n = 2a£ n /c and z = x n p. It is 
evident that x n = rn = 2irnT '/T e ff, where the effective 
temperature is defined by ksT e s = he/ (2a) In new 
variables Eq. (0) takes the form 
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where Qi,2 in (0) are expressed now in terms of x n , z. 

The sum in (g) can be calculated with the help of the 
Abel-Plana formula (l8|. The result is 
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The first term in the right-hand side of (Q) is the 
Casimir force at zero temperature, the second one takes 
into account the temperature corrections. The zero tem- 
perature contribution was calculated in [l9| numerically 
by the use of optical tabulated data for the complex re- 
fractive index (an alternative computation |pof contains 
some errors which are indicated in p9|). Independently, 
in it was determined by perturbation theory up to 
the fourth order in the small parameter 5q /a (Sq being 
the effective penetration depth of electromagnetic zero 
point oscillations into the metal). Thereby, the plasma 
model was used for the dielectric permittivity 



^0 = 1 + ^ = 1 + ^1 
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where u> p is the effective plasma frequency, and uj p = 
2auj p /c, so that a = 1/uj p = So /(2d). The results of 
|l9f and Q are in good agreement for space separations 
a > A p = 2nc/ujp. It is well known that the plasma 
model does not take into account the contribution of re- 
laxation processes which are taken into consideration by 
the Drude model (see below). However, the variation 
of the Casimir force obtained by both models remains 
smaller than 2% §§. 

Let us calculate the second term of (||) in the appli- 
cation range of plasma model and under the condition 
T <C T e g. To do this we use the representation of (||) 
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Introducing xq < 1 one can write 
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Considering firstly the case i — 2, we notice that in the 
plasma model (||) the second term from the right-hand 
side of (0) does not depend on x. Expanding Q^ 1 from 
the first term of (ffi) into a series in powers of z and inte- 
grating one arrives at the result 
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where C = const. 

The even powers of x evidently do not contribute to 
the second term of (^) . As a result there is only one tem- 
perature correction originating from Q2 which is caused 
by the term x 3 /6 and which does not depend on ui p . Sub- 
stituting (||) into the second term of @ one obtains 
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P p (a) = —Tr hc/(2A0a ). Note that we neglect the 
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corrections O [(T/T cfi ) 5 } . 

Consider now i = 1 in (Q). In this case both the first 
and the second terms in the right-hand side depend on 
x. The second term, however, is an even function of x 
and for that reason it does not contribute to (|]). Let us 
expand the quantity z 2 Q^/ 1 in powers of small parameters 
a and z. Integrating the obtained series between the 
limits z = x and z — xq < 1 we obtain 
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where the quantity 
not contribute to 
(T/T eff ) 5 or higher, 
term of (0) we get 

A T F$(a)=F$( a ) 



ifpj (x) contains terms which do 
) or lead to contributions of order 
Substituting (noh into the second 



T 

T e f£ 



30C(3) So (T_ 

. Teg 



(11) 

where £(3) w 1.202 is the Riemann zeta function. 

Now, let us take together ([)]), ( pi] ) and the zero tem- 
perature contribution given by the first term of (^). In 
pl| the last one was calculated up to the fourth order. 
Here we add two more orders. The final result is 
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where C2 = 24, and the other coefficients are 
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For 5q = (perfect conductor) Eq. ( |12[) turns into the 
well known result [13-15]. It is significant that the first 
correction of mixing finite conductivity and finite temper- 
ature is of order (T/T c g) 3 , and there are no temperature 
corrections up to (T/T Q g) 4 in the higher conductivity cor- 
rections from the second up to the six order. 

Analogous calculations can be performed for the con- 
figuration of a sphere (lens) of radius R above a plate 
starting from the force 
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This formula is obtained from (|]J) using the proximity 
force theorem ^2j, Q1.2 arc defined in (j|). After straight- 
forward calculations, using pl| ] for the zero temperature 
contribution, the result is 
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where i^ 0) (a) = -7r 3 fcJ?/(360a 3 ), c l = 3c*/(3 + i), c, 
are defined in ( |l3| ) . For the perfect conductor <5 — > the 
known asymptotic behavior |?J is reproduced. 

Now we consider space separations a for which T ~ T c g 
or even larger. In this case perturbation theory in T/T e g 
does not work. Let us compute the values of temperature 
force ( 14 ) numerically in dependence on a for Al surfaces 
used in experiments |^|| with lo p = 1.92 x 10 16 rad/s 
||, T = 300 K, and R = 100 fim. The numerical re- 
sults are shown in Fig. 1 by the solid curve. In the same 
figure the asymptotic behavior ( plj ) is presented by the 
pointed line. The dashed line shows the Casimir force 
at zero temperature (but with account of finite conduc- 
tivity). Here, the force was computed by Eq. ( |l4| ) in 
which the sum has been changed into the integral Jl5| . 
It is seen that perturbation theory works well within the 
range 0.1 fim < a < 3.5 fim (note that all six perturba- 
tion orders are essential near the left verge of this inter- 
val). Starting from a = 6 /im the solid line represents the 
asymptotics at large separations (temperatures) 
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This result follows from the term of ( |l4| ) with n = (the 
other terms being exponentially small in T/T e g). For 
So = one obtains from (jl^) the known expression for 
perfect conductors [7,13-15]. Finite conductivity correc- 
tions of higher orders do not contribute at large separa- 
tions. 

By way of example, consider the contribution of tem- 
perature correction for the Al sphere and the plate of 
experiments M M at smallest separations (a ~ 0.1 /im) as 
calculated in P%. According to the numerical results ob- 
tained by the plasma model AtF p i = F p i(T) — F p i(T = 
0) w 0.03 pN, where F p i (T) is computed by Eq. (||). 
Almost the same result is obtained by the perturbative 
result of Eq. ( )l5| ) . These values fall far short of the com- 
putational result \AtF p i\ »j 4pN presented in p7[ . Al- 
though the computational procedure is not described ex- 
plicitly in |17j we have been able to reproduce the value 
obtained there as follows. 

In jlTj , instead of the plasma model, the Drude model 
was used at small frequences for which the dielectric per- 
mittivity on the imaginary axis is 



e(iO - 1 + 



(17) 



Here, 7 is the relaxation frequency. Substituting the di- 
electric permittivity (17) into ( |l4| ) and performing com- 
putations the incorrect values \AtF p i\ « (2.5 — 8.5) pN 
are obtained when separations decrease from 0.1 /zm till 
0.09 fim (compare with « 4pN in jl?]]). It arises for the 
following reasons. 

When performing computations, proper allowance 
must be made for a critical issue resolved by J. Schwinger, 
L. L. DeRaad, Jr., and K. A. Milton |]|. In line with 
p5| the prescription should be adopted that we take the 
limit e — > 00 before setting £ = in order that the 
Casimir force between perfect conductors be obtained 
from Eqs. (|l]), (jlj). Otherwise the n — terms of (Q), 
(|l4|), containing Q2, would not contribute, which would 
imply incorrect limits both at low and high temperatures 
15 . For a real metal the prescription of [jl5| is satisfied 
automatically when theplasma model (||) is used. In the 
case of Drude model (|17|), however, the contribution of 
Q2 in the terms of (|T|), (|14[) with n = is identical zero. 
Because of this, it is impossible to follow the prescription 
of p5[ . It is necessary to stress that if the prescription 
of is not carried out for the calculation of the zeroth 
term of (|lj) one would obtain —((3)RkBT/(8a 2 ) instead 



of (16) (remind that the term containing Q2 does not 
contribute in this case). The last expression is evidently 
incorrect. It is two times smaller than the main con- 
tribution to (|l6|), which is valid for perfect conductor, 
and is independent of the actual value of conductivity 
of the real metal under consideration. Now it has been 
evident that the extra contribution of «4pN discussed 
above originates from the missing contribution of Q2, 
when the Drude model is used, and is equal to it by 
the modulus. Because of this, Drude model can not be 
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use d to provide a correct extraction of the case of metals [17] V. B. Sveto voy and M. V. Lokhanin, e-print quant 



from the Lifshits formula for dielectrics. 
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In conclusion it may be said that the joint effect of non- 
zero temperature and finite conductivity on the Casimir 
force was examined. It turned out to be in agreement 
with the previous knowledge for the real metals at zero 
temperature from one side and for the perfect conductors 
at non-zero temperature from the other. (Note that some 
of the above results related to the plasma model only were 
obtained independently in the recent preprint f24{| .) The 
obtained results are the topical ones for the interpretation 
of precision measurements of the Casimir force. 
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FIG. 1. The Casimir force in pN as a function of the 
surface separation in configuration of a sphere above a disk. 
The solid line represents the computational results obtained 
by Eqs. (|l3). The dotted line is calculated by the per- 
turbative Eq. (|15|) up to sixth order in relative penetration 
depth and fourth order in relative temperature. The dashed 
line is the zero temperature result. 
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